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Transformations of polynomial ensembles 

Arno B. J. Kuijlaars 


Dedicated to Ed Saff on the occasion of his 70th birthday 


Abstract. A polynomial ensemble is a probability density function for the 
position of n real particles of the form Ylj<k( x k ~ x o) det [fk( x j)Yjk— i> 
for certain functions /i,..., fn- Such ensembles appear frequently as the joint 
eigenvalue density of random matrices. We present a number of transforma¬ 
tions that preserve the structure of a polynomial ensemble. These transforma¬ 
tions include the restriction of a Hermitian matrix by removing one row and 
one column, a rank-one modification of a Hermitian matrix, and the exten¬ 
sion of a Hermitian matrix by adding an extra row and column with complex 
Gaussians. 


1. Polynomial ensembles 

A polynomial ensemble is a probability density function on K™ of the form 
(1-1) V{x!,. ..,x n ) = 2 - A„(x) det [fk{xj )] n jk=1 , 

where /i, ...,/„ is a given sequence of real-valued functions, 

A n (x) = - Xj) = det [:Tj _1 ]yfc=ij x = (aq,. ..,x n )e K”, 

j<k 

denotes the Vandermonde determinant, and Z n is a normalization constant. Cer¬ 
tain conditions on the functions fi ,..., /„ have to be satisfied to ensure that m 
is indeed a probability density. For example, the functions should be linearly inde¬ 
pendent and the integrals 


id - 1 


fk{x)da 


j,k = l,...,n 


should be convergent. In addition, m has to be non-negative for all possible 
choices oi (xi,... ,x n el". The probability density (11.11) only depends on the linear 
span of the functions f ±,..., f n , as one can see by applying column transformations 
to the second determinant in (ED- 

A special case arises if /*,( x) = x k 1 w(x), for k = 1,... n, with w an integrable 
non-negative function on R such that the moments up to order 2n — 2 exist. In 


l 
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that case ed can be written as 


(1.2) V{x-L,...,x n ) = — Al( x ) Y[w(xj), x = (x \,..., x n ) 

n 1=i 


and ed is known as an orthogonal polynomial ensemble HU, as the analysis of 
ED relies on the polynomials that are orthogonal with respect to w. The ensem¬ 
bles ed arise as the joint probability density of eigenvalues of unitary invariant 
ensembles of Hermitian random matrices [9]. The polynomial ensembles ED also 
include the multiple orthogonal polynomials ensembles, see m, where also more 
examples from random matrix theory are given. 

On the other hand, we have that ED is a special case of the more general class 
of biorthogonal ensembles, see im 


(1.3) 


V^xx ,... ,x n ) = — det [gk{xj)]l k=1 det [/fc(zj)]" )fc=1 , 


which involves two sequences gi,...,g n and /i,...,/ n of given functions. It is 
known that ED is determinantal [Tj, which means that there exists a kernel I\ n : 
MxIaM such that 

V , - - *, x n ) — det \I\. n {xj , £/c )]j f^—i 

and such that for every m = 1,..., n — 1, 

f P{xx ,..., 3?n)dx m _|_i • • • dx n — j - det [A. n (xj , £&)] ~ ^—x • 

Ir-» n - 

The correlation kernel K n has the form 

n 

K n (x,y) = y^ipk(x)(j)k(y) 

fc=l 

where spanj^i,..., x/j n } = spanj^i,..., g n }, span{^i,..., </>„} = span{/i, ...,/„} 
and 

/ OO 

ipj(x)(j)k{x)dx = Sj,k j,k = l,...,n. 

-OO 

In the case (ED we may write 

n— 1 

K n (x,y ) = y2, P oi x )Qo{v) 
j =0 

where Pj is a monic polynomial of degree j for j = 0,..., n — 1, the dual functions 
Q o,..., Qn-i are in the linear span of /i,..., /„, and the biorthogonality condition 


(1.4) 


/ OO 

Pj(x)Q k {x)dx = 5j,k, j, k = 0,... ,n - 1 

-OO 


is satisfied. In this case we can also consider the monic polynomial P n of degree n 
such that El also holds for j = n. This polynomial is given by 


t=i 


Pn{x ) = E 
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where the averaging is over {x\,...,x n ) in the polynomial ensemble (11.11) . If the 
points x \,..., x n come from eigenvalues of a random matrix, then P n is the average 
characteristic polynomial. 


2. Known transformations 

This paper discusses a number of transformations that preserve the structure of 
a polynomial ensemble. These transformations come from random matrix theory, 
and the typical setting is the following. We assume that X is a random matrix whose 
eigenvalues (or squared singular values) are distributed according to a polynomial 
ensemble (EH)- Then we perform a certain transformation to obtain from X a new 
random matrix Y. and the result is that the eigenvalues (or squared singular values) 
of Y are again a polynomial ensemble. 

2.1. Product with Ginibre matrix. The first example of such a transfor¬ 
mation comes from recent work of the author with Dries Stivigny 1201 . It deals 
with the squared singular values of rectangular matrices. Recall that the squared 
singular values of a rectangular complex matrix X are the eigenvalues of X*X. 
The transformation on X is multiplication by a complex Ginibre matrix, where a 
complex Ginibre matrix is a random matrix whose entries are independent standard 
complex Gaussians. 

Theorem 2.1. Let n, l , v be non-negative integers with 1 < n < l. Let G be an 
(n + u) x l complex Ginibre matrix, and let X be a random matrix of size l x n, 
independent of G, such that the squared singular values x \,..., x n are a polynomial 
ensemble (ED for certain functions fi,...,f n defined on [0, oo). Then the squared 
singular values y\,... ,y n ofY = GX are a polynomial ensemble 


(2.1) 

n(y ) det [5fc(%)]", fe= 1 , 

^n 

all yj > 0, 


where 




(2.2) 

9k{y) = [ x' J e~ x f k —. 

Jo \xJ x 

y> o. 


Proof. 

See [20], where the proof is based on 

ideas taken from [3] [4]. 

□ 


Note that gu in (|2.2I) is the Mellin convolution of cc i—>- x v e x with fk- 
Theorem EH can be applied repeatedly and it follows that the multiplication 
with any number of complex Ginibre matrices preserves the structure of a polyno¬ 
mial ensemble for the squared singular values. 

Theorem l2.1l was inspired by earlier results by Akemann et al. ia 1 1 on products 
of random matrices. In these papers the authors considered products of complex 
Ginibre matrices (that is, X is also a complex Ginibre matrix) and they obtained 
the structure ED ED. where in this case the functions gk in ED are expressed 
as Meijer G-functions. This result has since then been used in [211(22] to determine 
the large n scaling limit of the correlation kernel at the hard edge, and in [2] to 
calculate the Lyaponov exponents as the number of matrices in the product tends to 
infinity. See also (Hum in] for other recent results on singular values of products 
of random matrices. 


4 


ARNO B. J. KUIJLAARS 


2.2. Product with a truncated unitary matrix. Theorem 12. ll has an ex¬ 
tension to a product with a truncated unitary matrix. A truncation T of a matrix 
U is a principal submatrix of U . We assume that U is a Haar distributed random 
unitary matrix and then T is also a random matrix. 

Theorem 2.2. Let n,m,l,v be non-negative integers with n < l < m and 
m > n + u + 1. Let T be an (n + v) x l truncation of a Haar distributed unitary 
matrix U of size m x m. Let X be a random matrix of size l x n, independent of 
U, such that the squared singular values Xi,..., x n of X are a polynomial ensemble 
(ED for certain functions fi, ■ ■ ■, f n defined on [0, oo). Then the squared singular 
values y\,... ,y n ofY = TX are a polynomial ensemble 

(2-3) -i- a n (y) det [3fc(2/j)]" |fc=1 , all yj > 0, 

where 

(2.4) 9 k(y)= f X V (1- X ) m - n - v - 1 fk (?-)—, y> 0. 

Jo X 

Proof. See m- □ 

If we let ui —y oo in Theorem l2.2l then \JrnT tends in distribution to a complex 
Ginibre matrix. Also (1 — —tends to e~ x as m —> oo. In this way 
Theorem m can be obtained as a limiting case of Theorem 12.21 

Theorems o and 12.21 can be used repeatedly and it follows that the squared 
singular values of a product of any number of Ginibre matrices with any number 
of truncated unitary matrices are a polynomial ensemble. 

2.3. Overview of the rest of the paper. Inspired by these results we give 
an overview of other transformations that preserve polynomial ensembles. The 
transformations are based on known random matrix theory calculations, see m 
[15] . and our aim here is to emphasize the interpretation as a transformation of 
polynomial ensembles. 

The first such transformation comes from matrix restrictions. Here we are 
working with a Hermitian matrix A' and we remove one row and one column to 
obtain Y. If X is random with eigenvalues that are distributed as a polynomial 
ensemble then the eigenvalues of Y are also distributed as a polynomial ensemble. 
This is our first result, see Theorem 13.21 The proof relies on a fundamental result 
of Baryshnikov [ 6 ] , see Theorem 13.11 below. 

Then we extend this to the situation where A is a positive semidefinite matrix 
with a fixed number of zero eigenvalues. Again we find that matrix restriction 
for random matrices of this type leads to a transformation result for polynomial 
ensembles, see Theorem 14.21 Interestingly enough, we can make a connection with 
the product with a truncated unitary matrix, as we find in this way an alternative 
proof for Theorem 12.21 

In Section [5] we consider a transformation from X to Y = X + vv* where X 
is Hermitian, and v is a column vector of independent complex Gaussian entries. 
This rank-one modification is also a transformation of polynomial ensembles as we 
show in Proposition 15.II The argument is based on a result of m- 

Finally, in Section [6] we consider a transformation where we extend the Hermit¬ 
ian matrix X by adding an extra column v with independent complex Gaussians, 
and an extra row (y* c) consisting of v* and a real number c that has a real 
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normal distribution. Under appropriate conditions on the variances, we again find 
a transformation of polynomial ensembles, see Proposition 16.21 This is based on 

mini- 


3. Matrix restrictions 


Let X be annxn Hermitian matrix with distinct eigenvalues x\ < X 2 < ■ • • < 
x n . Let U be a Haar distributed unitary matrix of size n x n and let Y be the 
(n—1) x (n—1) principal submatrix of UXU* with eigenvalues y\ < y -2 < ••• < y n - 1 - 
With probability one we have strict interlacing of eigenvalues 


(3.1) xi < yi < x 2 < yi < • • • < y n -i < x n . 

The following theorem is due to Baryshnikov (reformulation of |6j Proposition 4.2]). 

Theorem 3.1. If X and, Y are as above, then the (random) eigenvalues yi,..., y n -i 
of Y have the joint density 


(3.2) 

on the subset o/M” _1 


(n — 1)! 


A„i(y) 


A n (i) 

defined by the inequalities (EH). 


The interlacing condition is expressed by the determinant 


(3.3) 


det \Xx k <yj\ j k=l > 


1 if x < y 
0 otherwise, 


with y n := +oo. Indeed, for all mutually distinct values x k and yj, the determinant 
in EH is 1 if and only if the interlacing condition holds and it is zero otherwise. 
The determinant in (13.31) has all ones in the last row. We can reduce it to an (n— 1) x 
(n — 1) determinant by subtracting the last column from every other column, and 
expanding along the last row. This results in the determinant det [x*fc<j/^<®n]"fc= 1 - 
It means that the density (13.21) can be written as 


(3.4) 


An-i(y) 

A„(x) 


det [Xx h < VJ ]l k=1 


An-i(y) 
A „(x) 


det \xxk<yj<x„ 


n— 1 
j,k= 1 ’ 


y = (yi, ■ ■ ■ ,2/n-i) e R 71 1 , y n -=+oo, 


where (13.41) is considered as a probability density on the unordered eigenvalues of 
Y, i.e., as a probability density on R" -1 . This accounts for the disappearance of the 
factor (n — 1)! from (13.21) . Note that (13.41) is a polynomial ensemble with functions 
y Xx k < y <x n , for k = 1,2, 1. 

Let us now assume that X is random, independent of U, and that the eigenval¬ 
ues of X are a polynomial ensemble. Then the eigenvalues of Y are again a poly¬ 
nomial ensemble. For this it is important that the normalization constant A ^ in 
(13.41) depends on X via the Vandermonde determinant A n (x) in the denominator. 
We also need the Andreief identity, see uni Chapter 3], 


(3.5) J det[$ k (xj)]™ >k=1 det[ip k (xj)]™ k=1 dy,(xi) ■ ■ ■ dy(x n ) 

4>j(x)ip k {x)dn(x) 


= n\ det 


j,k =l 
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where y is a measure on a space X, n £ N, and (f>i, ..., (j> n , tfi, ... , tj) n are arbitrary 
functions such that the integrals converge. The result is the following. 


Theorem 3.2. Suppose that X is a random n x n Hermitian matrix whose 
eigenvalues are a polynomial ensemble (ED with certain functions /i,..., f n . Let 

Y be the principal submatrix of UXU* of size (n — 1) x (n — 1), where U is a Haar 
distributed unitary matrix, independent of X. Then the eigenvalues yi,... ,y n -1 of 

Y are a polynomial ensemble 

(3-6) -J— A„_i (y) det [gk(yj)]jk=i, 

Z n -1 

with 

(3.7) 9k(y)=f fk{x)dx , k = l,...,n-l, 

J — oo 

where /i,..., f n -i are a basis for the vector space 

n pOO 

(3.8) {/ = ^ Cjfj | ci,,.., c„ e M, / f(x)dx = 0}. 

i=i J ~°° 

Proof. From (IQ) it follows after averaging over the polynomial ensemble 
ED that the eigenvalues of Y have joint density 

{y) J det [xx k < yj ]l k=1 det [/fc(^)]" fc=1 dx i • • • dx n 

with y n := +oo. Because of (13.511 we find that this is 


(3.9) 


-A„_i (y) det 


rvj 


fk(x)dx 


j,k= i 


Since y n = +oo, the last row of the second determinant in m contains the 
constants fk(x)dx. Change from fi,...,f n to another basis fi, ■ ■ ■ f n where 

/i,..., f n -1 belong to the subspace (13.811 . Then after performing suitable column 
operations we obtain the density 


Z n -\ 


A n _i (y) det 


rvj 


fk{x)dz 


j,k =1 


Then f_ f n (x)dx ^ 0, since otherwise the full last row in the determinant would 
be zero. By expanding the determinant along the last row we obtain (13.611 with 
functions (ED and a possibly different constant Z n -\. □ 


An analogous result holds for singular values. 


Corollary 3.3. Suppose that X is random (n + v)x n matrix whose squared 
singular values are a polynomial ensemble ED- Let Y be the (n + u) x (n — 1) 
principal submatrix of XU where U is Haar distributed unitary matrix, independent 
of X. Then the squared singular values yi,..., y n -\ ofY are a polynomial ensemble 

(3.10) -J— A„_i (y) det [gk(yj)]^t i , all Vj > 0, 

^n -1 

with 

(3.11) gk(y)=[ fk(x)dx , fe = 1,..., n — 1 

Jo 
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where f 1 ,.. •, f n -1 are a basis for the vector space 
n /»00 

(3.12) {/ = I ci,...,Cn € K, / /(ai)da; = 0}. 

7=1 

PROOF. We can apply Theorem 13.21 since Y*Y is the principal submatrix of 
size (n — 1) x (n— 1) of X*X. The integration in (13.111) and (13.121) starts at 0 since 
the functions are defined for x > 0 only. □ 


4. Restrictions of positive semidefinite matrices 

The following is a variation on Theorem l3.ll It can also be obtained as a special 
case of m Corollary 1]. 


Proposition 4.1. Let m > n + 1 and let X be an m x m positive semidefinite 
Hermitian matrix with n simple non-zero eigenvalues 0 < X\ < X 2 < • • • < x n and 
an eigenvalue 0 of multiplicity m— n > 1. Let Y be the (to — 1) x (to — 1) principal 
submatrix of UXU* where U is a Haar distributed unitary matrix of size m x to. 
Then with probability one, Y has exactly n non-zero eigenvalues 0 < y\ < 1/2 < 
• • • < y n that satisfy the inequalities 

(4.1) 0 < yi < xi < y 2 < x 2 < ■ ■ ■ < x n -± < y n < x n , 

and these non-zero eigenvalues have the joint density 

u 2 ) (™-i)! n yr n - 1 M) 

(to — n — 1)! x m-n a„(x) 

restricted to the subset of y = (y ±,..., y n ) £ K" defined by the inequalities (14.11) . 


PROOF. For m = n + 1 this follows immediately from Theorem 13.11 and so 
we assume in the proof that to > n + 2. We approximate X by a matrix A with 
eigenvalues 01 < ■ • • < a m - n < x\ < ■ ■ ■ < x n with aj’s close to zero. Let B be the 
principal submatrix of UAU* of size (to — 1) x (to — 1), which with probability one 
has distinct eigenvalues 61 < • • • < b m -n- 1 < yi < • • • < y n that interlace with the 
eigenvalues of A. By Theorem 13.11 the joint density of these eigenvalues is 


(4.3) 


^ | A m -n-l(b) rij-,fc(3/fc bj) A n (y) 

A m _„(a) aj) A n (x) 


on the subset of R m given by the interlacing relations. The induced density on 

yi,---,y n is 


(4.4) (to 


1 )! 




dbm.— n— 1 


A m —n— 1 ( 6 ) Ylj,k(yk bj) 
A m-n{a) rij,fc( Xfe — a i) 


A n(y) 
A n{x) 


In the limit where all aj —> 0 , j = l,...,m — n, we also have bj —> 0, j = 
1,... ,m — n — 1. Then the factors ■ k {y^ — bj) and k {xj- — aj ) in (14.41) tend to 
IL V k~ n ~ l and n fe x™ ", respectively. The resulting m — n — 1 fold integral can 
be evaluated as 


(4.5) 




A m _ n _i(5) 
A m—n[a) 


1 

(to — n— 1 )! 


The result is the joint density (14.21) for 

□ 


and this does not depend on 01 ,... a m - n 
the non-zero eigenvalues of Y. 
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The inequalities m are encoded by the determinant 

det [xo <yj<x k \j k=1 

which for strictly increasing yi < 2/2 < • • • < y n is 1 if the interlacing ED holds 
and 0 otherwise. Then (Id.21) can be alternatively written as 

(to — 1)! ” y™-"- 1 A n (y) 


(4.6) 


i!(m — n — 1 )! 


n 

k= 1 


A n (x) 


det [x. 


0< yj <x k \ j k=1 


which is now considered as a density on [0,oo) n for unordered eigenvalues. Note 
that m is a polynomial ensemble on [0, 00 ) with functions y i-A y m n 1 Xo<y<x k 
for k = 1 ,..., n. 


Theorem 4.2. Let n < to—1 and v < to — n — 1 be positive integers. Let X be a 
random positive semidefinite Hermitian matrix of size mxm with a zero eigenvalue 
of multiplicity m — n > 1 and non-zero eigenvalues x\ ,..., x n that are a polynomial 
ensemble ED for certain functions fi,..., f n on [0, 00 ). Let Y be the principal 
submatrix ofUXU* of size (n + u) x (n + v), where U is a Haar distributed unitary 
matrix, independent of X. Then, with probability one, Y has exactly n non-zero 
eigenvalues yi, ... ,y n , and these non-zero eigenvalues are a polynomial ensemble 

(4-7) -=~ A n(y) det [gk(yj)]" k=1 all yj > 0, 

^n 

where 

( 4 - 8 ) g k (y) = f x v (l - x) m - n ~ v - :1 f k (?) —. 

J 0 x 

Proof. We first assume that v = m — n — 1. Then Y is obtained from UXU* 
by removing one row and column and we can apply Proposition l4.fl and in particular 
its reformulation in (14.61) . Averaging (14.61) over the polynomial ensemble (11.11) we 
obtain the joint density 

1 f tu w m_n_1 „ 

-jpA n (y) / Yl m-n det [A-(^)]” fc=1 det [xo< w <x*]" fc=1 dx 1 ■ ■ ■ dx n , 

^n J[0,oo)” fc—'y x k 


for a certain constant Z n (which also depends on in). By the Andreief identity 
ED this leads to 


-~~A n (y) det 
Zn 


{•00 m—n—1 

a 


—fk{x)xo< Vj <xdx 


with a new constant Z n . This is a polynomial ensemble (14.71) with functions 

POO — TL — 1 

gk(y) = / m-n fk( x )X0<y<x dx 

J 0 X 

/ oo _ .m—n—1 

x m- n Mx)dx, y> 0. 

The substitution x - in (14.91) leads to the expression (14.81) with v = to — n — 1. 
This is the Mcllin convolution of f k with the function Xm-n-i where we define 

Xk : [0, 00 ) —>■ R : x x k Xo<x<i ■ 
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Thus gk = Xm-n-i * fk, if v = m — n — 1, where * is used here for the Mellin 
convolution 

( f*g)(y)=[ 

Jo ^ x/ x 

For general v < m — n — 1 we can use the above argument repeatedly, and 
we find a polynomial ensemble m with functions gk that are iterated Mellin 
convolutions of the functions fk, namely 

9k = Xv * Xv+1 *■■■* Xm—n—l * fk, k = 1, . . . , Ti¬ 
lt is easy to calculate that 

(Xv * Xv+1 *■■■* Xm—n—l') {x) = x"{l - x) m ~ n ~ V ~ 1 , 0 < X < 1, 

and thus we obtain the formula (14.81) for the functions gk- □ 

An attentive reader may have noticed that the formula for gk in (14.81) coincides 
with the one appearing in (12.41) in Theorem 12.21 This is no coincidence since we 
can use Theorem 14.21 to give an alternative proof of Theorem 12.21 

Proof of Theorem 12.21 Let X be an l x n matrix, and put 

^ _ (XX* 0 
0 0 

which is an m x m matrix with m — l rows and columns containing only zeros. It is 
clear that the squared singular values of X are equal to non-zero eigenvalues of X. 

Also if U is a unitary matrix of size m x m and T is its left upper block of size 
(n + v) x l then 

TXX*T* *\ 

* *1 

where * denotes a certain unspecified entry, whose value is not important for us. 
In other words, TXX*T* is equal to the principal submatrix of UXU* of size 
{n + v) x (n + u). Assuming U is Haar distributed over the unitary group, and 
the squared singular values of X are a polynomial ensemble (ED, independent 
of U , we then find from Theorem 14.21 that the non-zero eigenvalues of TXX*T* 
are a polynomial ensemble ED with functions (14.81) . The non-zero eigenvalues 
of TXX*T* are the same as the squared singular values of TX and Theorem 12.21 
follows. □ 


UXU* = 


5. Rank one modification 


Let X be a Hermitian n x n matrix with eigenvalues x\ < X 2 < • • • < x n . 
We take Y = X + vv* where v is a vector of length n. Then the eigenvalues 
Uj of Y interlace with those of X , as follows from the Courant-Fisclier Theorem, 
see e.g. }23l chapter 7.5]. We let v = (vi,...,v n y be a vector of independent 
complex random variables whose real and imaginary parts are independent and 
have a AT(0,1/2) distribution. Then the distribution of the eigenvalues of Y is 
given in m Appendix E] as 


(5.1) 


3 = 1 


A n (x) 


on the subset of R™ given by the interlacing conditions 


( 5 . 2 ) 


xi < 2/1 < x 2 < ■ ■ ■ < x n < y n 
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The following result is an immediate consequence. 


Proposition 5.1. Let Reu,-, Imu,-, for j = 1,... ,n be mutually independent 
normal random variables with mean zero and variance 1/2. Let X be a random 
Hermitian matrix of size n x n, independent of v = (iq,... whose eigenval¬ 
ues are a polynomial ensemble JED with certain functions /i . Then the 

eigenvalues yi,... ,y n ofY = X + vv* are a polynomial ensemble 

(5-3) J- A n (y) det \g k (.yj )}] fc=1 , 

Z n 

where 

POO 

(5.4) g k ( y )= e~ x f k (y - x)dx, fc = l,...,ra. 

Jo 

Thus gk is the convolution of //- with x K > e~ x \x>o- 


PROOF. The interlacing (15.21) is encoded by a determinant, and it follows that 
<E3D is a polynomial ensemble 


(5.5) 


1 " 

sir 

i=i 


-(.Vj-xj) 


An(y) 

A„(ar) 


det [y : 


x k<Vj J j } /c=l ’ 


where now we disregard the ordering of the yj ’s and consider (15.51) as a probability 
density on K™. 

We average over X\,..., x n distributed as in JED- By Andreief’s identity (13.51) . 
we obtain for the density of the eigenvalues of Y 


■=-A n (y) det 

Changing variables x yj — x in the integral in the determinant, we arrive at (15.31) 
with functions (15.41) . □ 


rvj 


,-{Vj-x) 


fk{x)dx 


J j , k =1 


Here is a variation on the same theme. 


Proposition 5.2. Let n,v > 1. Let Reu^, lm . Vj , for j = 1 ,...,n + u be 
mutually independent normal random variables with mean zero and variance 1 / 2 . 
Let X be a random (n + v) x (n + u) positive semidefinite Hermitian matrix, in¬ 
dependent of v i,..., v n , with exactly n positive eigenvalues nq,..., x n that are a 
polynomial ensemble (0) with certain functions fi,-..,f n on [0,oo). Then, al¬ 
most surely, Y = X + vv* has an eigenvalue zero of multiplicity v — 1 and n + 1 
positive eigenvalues y i,... ,y n + 1 that are a polynomial ensemble 


(5.6) 


-J— A n+1 (y) det [gk(yj)]™th > 

^n+1 


with functions 

9i(y) = y’'~ 1 e ~ v , 

(5.7) . _ fy 

gk+i(y) = y v e v X v e x fk{x)dx , 

J c 


all yj > 0 , 


k = 1 ,... ,n, 


for y > 0 , where c € ( 0 , oo) is an arbitrary but fixed positive real number. 


We may also take c = 0 or c = oo in provided that the integrals are all 
convergent. 
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PROOF. We approximate X by A with distinct eigenvalues ai < ■ ■ ■ < a u < 
x\ < ■ ■ ■ x n where the a,j are close to 0. Then B = A + vv* has eigenvalues 
bi < ■ ■ ■ < &j,_i < yi < ■ ■ ■ < y n + 1 that interlace with those of A , with a joint 
density, see ED, 


n e j li e ~ yj li f "' ii < :j: ‘ ——- l — 

3 = 1 3 = 1 t =1 t=l 


+«) (n-=i nLi(*fc - a t)) A »( a; ) 


subject to the interlacing conditions. 

We restrict this to the y- variables by integrating out bi, ..., 1 - This gives 

the joint density for y u ..., y n+1 


n+1 n v 


li n n 


_ A »+i (y) _ 

(rij=l Ylk=l( x k - %)) A n(z) 




v-1 

db v -i ]^[ e~ bj A„_i(6) 

3 =1 


' v— 1 n+1 

n ~ 

V 1 = 1 fc=1 


In the limit where all aj ->0we also have that all bj —> 0 because of the interlacing. 
Then A — > X 1 B — > Y, and using also (14.51) we find the limiting joint density for 
the nonzero eigenvalues yi,, y n +i of Y 


n+1 


-'n+l 


n (#r'< 


-Vj 


3 =1 


)n 

3 =1 


(z-^) 


A n +1 (y) 

A n (x) 


subject to the interlacing 0 < y\ < x\ < ■ ■ ■ < x n < y n + 1 - The interlacing is 
encoded by the determinant det [Xyj<x k <y n+1 ] 

Next, averaging over the polynomial ensemble m and using the Andreief 
identity m , we find in a now familiar fashion a joint density 


(5.8) 




A n+1 


(y) det 


n+1 


y--^~ Vj 


x v e x fk{x)x 


Vj<X<Vn +1 


dx 


J j,k =1 


for the eigenvalues of Y. The n x n determinant in (15.81) is extended to an [n + 
1) x (n + 1) determinant by adding first a row with zeros and then a column with 
ones. Then after elementary column operations we easily arrive at the polynomial 
ensemble (15.61) with functions (15.71) and a possibly different constant Z n+ \. □ 


The two Propositions 15.11 and 15.21 have the following consequences regarding 
squared singular values of an extension of a matrix by one row or one column. 


Corollary 5.3. Suppose v > 0. Let X be an (n + u) x n random matrix with 
squared singular values 0 < X\ < • • • < x n that form a polynomial ensemble (ED 

with certain functions /i,. .., f n on [0, oo). Let Y = ( * J with v a random vector 

of independent complex Gaussians as in Proposition 1 5. 11 which is independent of 


X. Then the squared singular values y i, 
m with functions 


..., y n of Y are a polynomial ensemble 


9k{y)= [ e x f k (y - x)dx. 
Jo 


(5.9) 
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PROOF. The squared singular values of A' are the eigenvalues of A*A. The 

squared singular values of Y are the eigenvalues of (A* v) = A*A + vv*. 

Thus the result follows from Proposition 15. II The integration in (15.911 extends to y 
only, and not to oo as in (15.411 , since fk(x) is defined for x > 0 only, and we consider 
fk{y — x) to be zero if x > y. 

□ 

Corollary 5.4. Suppose v > 1. Let X be an (n + v) x n matrix with squared 
singular values 0 < X\ < ■ ■ ■ < x n that are a polynomial ensemble (ED with 
certain functions fi, ■ ■ ■, f n on [0, oo). Let Y = (A v ) with v a random vector of 
independent complex Gaussians as in Proposition 1 5. HI which is independent of X. 
Then the squared singular values y i,..., y n +1 ofY are a polynomial ensemble (15.611 
with functions (USD- 

Proof. The squared singular values of A are the non-zero eigenvalues of AA*, 
and the squared singular values of Y are the non-zero eigenvalues of YY* = AA* + 
vv*. Thus the result follows from Proposition 15.21 □ 

It is interesting to note that a combination of Corollaries 15.31 and 15.41 leads to 
the proof of one of the classical results of random matrix theory mm , namely 
that the squared singular values of a complex Ginibre matrix are distributed as a 
Laguerre ensemble. See also m Chapter 4.3.3] for a similar approach, and [12: 
for related results. 

Corollary 5.5. Suppose X is an m x n random matrix such that ReA^j, 
Im Ajj, i = 1,..., m, j = 1,..., n are independent normal random variables with 
mean zero and variance 1/2. Suppose v = m — n > 0. Then the squared singular 
values Xi ,..., x n of X have the joint density 

1 " 

(5.10) — ^(ar,- — Xi) 2 x^e~ Xi , all Xj > 0. 

n i<3 t=i 

PROOF. We use induction. It is easy to check Corollary [53] for m = n= 1. 

Assume Corollary 15.51 holds for certain m,n > 1. Note that (15.101) is a poly¬ 
nomial ensemble with functions fk(x) = x v+k ~ l e~ x for k = 1 Then by 

Corollary 15.31 it will follow that Corollary 15.51 also holds for m + 1 and n, and by 
Corollary [53] it holds for m and n + 1, provided that m > n. The calculations are 
straightforward and we do not give them explictly here. □ 


6. Matrix extensions 


In this final section we start from an n x n Hermitian matrix A and we are 
going to extend it to an (n+1) x (n+1) matrix by adding one row and one column. 
We write 


( 6 . 1 ) 



where v = (v\ V 2 • • • is a complex column vector, and c€ R is real. 

The following result was given by Forrester [12] and Adler, Van Moerbeke and 
Wang [1], see also ITT Chapter 4.3.2] and [16 . section 3.1], where the focus is on 
the situation where A' is an n x n GUE matrix. 
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Theorem 6.1. Suppose c, Reu,-, and Im Vj for j = 1 are independent 

normal random variables with mean zero, where c has variance 1 and Re Vj, Im Vj 
have variance 1/2. Assume X has simple eigenvalues X\ < X2 < ■ • • < x n . Then 
with probability one, the ordered eigenvalues t/i < 1/2 < • • • < y n + 1 of Y are simple, 
and strictly interlace with those of X: 

(6.2) yi < xi < y- 2 < x 2 < ■ ■ ■ < y n < x n < y n+ \. 


In addition, the eigenvalues of Y have the probability density 

n+1 


(6.3) 


1 


V2n 


iR^n 


1,? A n+1 jy) 

A n(x) 


e 2 i 


1=1 1=1 

on the subset o/R " +1 defined by the interlacing condition & 


PROOF. See Lemma 1 in [TJ or Proposition 6 in |12l . □ 

As before, there is an immediate consequence of Theorem 16.11 to polynomial 
ensembles. 


Proposition 6.2. Suppose c, Ren,-, and Inin, for j = 1 are mutually 

independent normal random variables with mean zero, where c has variance 1 and 
Re^j, Im iij have variance 1/2. Suppose that X is a random Hermitian matrix of 
size n x n, independent of c and v, whose eigenvalues are a polynomial ensemble 
m with certain functions / 1 ,..., f n . Then the eigenvalues y 1 ,..., y n +1 ofY given 
by (ICTTD are a polynomial ensemble 

(6-4) -J— A n +i{y) det \g k {yj)]™Xi 1 > 

An +1 


with functions 


(6.5) 


9 i{y) = e " y2 , 

gk+i{y) = e~i y2 f ei x2 f k (x)dx, 


k = 1 


Jo 

PROOF. The eigenvalues of X are distinct with probability one. We order 
them, say x\ < X 2 < • • ■ < x n . 

We use an interlacing determinant as in (13.311 to write the density (16.31) as 


Zn-\- 1 1 1 1 1 

1=1 1=1 

with x n +i = +00 and a certain constant Z n+ 1 


II 

1=1 


which 


n+1 

i,fc=1 


1 


Z, 


n+1 

n 


e-^n eK 


(y) 


A, 


1 Vi/y J . r -\n 

det \Xvi<x h <v n+ x\ jik=sl • 


( 6 . 6 ) 

/ y I 1 

1=1 1=1 

Then averaging (16.61) with respect to the polynomial ensemble (11.11) and using the 
Andreief identity (13.51) we obtain for the density of j/i,..., y n + 1 , 


- ffcT jl 

-II e“^ J 2 A„+i(?/)det 

j—i 

We extend the last determinant to an (n + 1) x (n + 1) determinant by adding 
first a row with entries f y +1 ^ x fk(x)dx for k = 1 , ...,n and then a column 


fUn + l 


'Vj 


fk{x)da 
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(0 • • • 0 l)* of length n + 1. Then we subtract the last row from each of the 

other rows and change the sign in each of the entries in rows 1 up to n. Then the 
result is (16.41) with functions (16.51) and a possibly different constant Z n +\. □ 


As a special case, we consider the polynomial ensemble (ED with functions 


fk(x) = x 
(6.7) 


k—lg — Tjx 


for k = 1,..., n. This is the same as 


1 

Z n 


A„(s) 2 R ( 

i=i 


Then by (16.51) we get g\(y) = e 2 s' 2 and 

9k+i(y) = e - ^ 2 J x k ~ l dx = ^y k e~i y2 for k = 1 ,..., n. 

The prefactor L is immaterial and it follows from Proposition 16.21 that the density 
function for the eigenvalues of Y is 

n +1 

-- A n+ i(y) 2 R e~* y i. 

Z n + 1 

It is well-known that (IQ) is the density of eigenvalues of GUE random matrix 
S3 HI] and we conclude, as already noted in ns Chapter 4.3.2], that we can use 
Proposition 16.21 to give an inductive proof of this basic result of random matrix 
theory. 
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